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Abstract 

The flow of a liquid crystal around a particle does not only depend on its 
shape and the viscosity coefficients but also on the direction of the molecules. 
We studied the resulting drag force on a sphere moving in a nematic liquid 
crystal (MBBA) in a low Reynold's number approach for a fixed director field 
(low Ericksen number regime) using the computational artificial compress- 
ibility method. Taking the necessary disclination loop around the sphere into 
account, the value of the drag force anisotropy (F±/Fn = 1.50) for an exactly 
computed field is in good agreement with experiments (~ 1.5) done by con- 
ductivity diffusion measurements. We also present data for weak anchoring 
of the molecules on the particle surface and of trial fields, which show to be 
sufficiently good for most applications. Furthermore, the behaviour of the 
friction close to the transition point nematic<->isotropic and for a rod-like and 
a disc- like liquid crystal will be given. 

PACS: 61.30.Jf; 61.30.Cz; 83.85.Pt 
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I. INTRODUCTION 



Most of the applications of liquid crystals are connected to their flow properties, either 
while processed or in the application itself. The reorientation of the director field for example, 
which is used in electro-optical devices, is linked to internal flow. The fastest response times 
needed for a further development are limited by the friction effects. But although the basics 
of the hydrodynamics of liquid crystals were laid about thirty years ago, most of the problems 
connected with flow are still unsolved. This is mainly due to the anisotropy of the system 
and the non-trivial connection of the direction of the molecules and the velocity. 

A deeper insight in the hydrodynamics of liquid crystals and the connection between 
macroscopic and microscopic properties would allow to predict the behaviour of particular 
materials and therefore to design special liquid crystals to obtain certain characteristics 
required. 

Precision experiments are often difficult to perform, since many of the standard tech- 
niques do not work for these materials. It would be useful to have more independent methods 
of measuring the viscosity, than only a traditional shear flow. A further technique, the falling 
ball experiment, was solved for an isotropic liquid by Stokes. It consists of a ball falling 
down in a cylinder driven by the gravitational force and measuring its equilibrium velocity. 
The viscosity rj can then be determined by the well known Stokes' formula F D = —6nrriv, 
which gives the relation between the friction drag Fd, the radius r of the ball and its velocity 
v. 

For liquid crystals this problem gets another dimension since the drag force also depends 
on the geometry of the system. It is obvious that the drag force on the sphere is different 
for the two particular cases of flow and director parallel and flow and director perpendicular 
to each other. Using a liquid crystal of rod-like molecules it becomes clear that it is easier 
to move the particle parallel to the general director field, i.e. along the long axis of the 
molecules, than to move it perpendicular to the director, i.e. perpendicular to the long axis. 
In the general case (arbitrary angle between flow and director) this results in the fact, that 
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the drag force is no longer parallel to the line of motion. There is a further component 
perpendicular to it, the so called lift force, which moves the particle sideways (see fig.f]). 
It is worth mentioning that this force does not contribute to the dissipative losses in the 
system, an effect well known from other areas of physics like electro-dynamics (a charged 
particle in a magnetic field is forced to change its direction without loosing or gaining any 
kinetic energy). 

A further problem is the influence of the director field n(r) on the flow since it does 
not only change but it also gives a contribution to the dissipative losses in the systems. 
In particular, regions of high gradients of the director field result in higher resistance to 
the flow. Such regions are mainly found around disclinations, which are often unavoidable 
due to the geometry of the system. If we consider, for example, perpendicular boundary 
conditions on the surface of the sphere and a uniform director field far away from it, there is 
a disclination loop around the sphere (see fig.0), which is unavoidable for topological reasons: 
The surface of the sphere corresponds to a s — 1 point defect and since the overall defect 
of the system must be zero (the director field is uniform far from the particle) this defect 
must be balanced by the disclination loop. The energy of the ring is roughly proportional to 
its length, therefore it is favourable to have it as small as possible. On the other hand, the 
rigid boundary conditions at the surface of the particle push the ring away from the sphere, 
so that the final position is given by the balance of the two effects. 

The drag force is sensitive to the radius of this loop. In both limiting cases (flow and 
director parallel/perpendicular to each other) the resistance is increased, but the magnitude 
of the influence is quite different. For the director perpendicular to the velocity the flow 
is parallel to the ring and it acts like a plate moved in the liquid crystal. For velocity and 
(general) director parallel to each other the flow is perpendicular to the ring, which does 
not only increase the cross-section of high director gradients around the ring and the liquid 
crystal flow but it also has a further effect: a certain amount of the liquid crystal has to 
flow through the gap between ring and sphere, where the director lies in the plane of the 
ring and is therefore locally perpendicular to the direction of the flow. As a consequence, 
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the anisotropy of the system, i.e. the ratio of the drag forces, decreases with an increasing 
radius of the disclination loop. That means that stronger boundary conditions lower the 
anisotropy. 

The theoretical problem of a liquid crystal flowing around a body has been addressed 
before. Diogo ]]J assumed the velocity field around the sphere to be the same as for an 
isotropic fluid and calculated the drag force for different angles between the director and 
the velocity. Roman and Terentjev obtained an analytic solution for the flow velocity for 
a fixed uniform director field, by an expansion in the anisotropy of the viscosity. Recently 
a group around Kneppe and Schneider gave solutions for the velocity of the liquid crystal, 
assuming a uniform director field, independent of the flow ||. 

All these solutions have their deficiencies. None of them, for instance, considered the 
distribution of the director field due to the boundary conditions on the particle. This will 
be done in this article, where the results are also compared with various approximations for 
the director n(r). 

The article is organised as follows: After a brief introduction to the basic equations of 
the hydrodynamics of liquid crystals which are needed in the next section, we give a short 
description of the numerical method we used to solve the equations of motion. Section IV 
gives the director fields we used and explains the limits when they are valid. The results for 
the drag properties and a comparison with experimental data are given in section V and, 
finally, we conclude with a discussion of possible experiments are given in chapter six. 

II. BASIC CONCEPTS 

In this section we give a brief summary of the nematic hydrodynamics that are used in 
this work. For derivations of these equations we refer the reader to the basic textbook [|J, 
see also ||. The stress tensor of a nematic liquid crystal consists of three contributions. 
They are the hydrodynamic pressure p, the viscous stress given by the tensor 

a' i:j = axniTijUkiiiAki + a 2 n j N i + a^Nj + a 4 A i:j + a 5 njn k A ik + a 6 nin k A jk . (1) 
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(here and below in this article we use the tensor index notation, i.e. an index appearing twice 
in a product means a summation over this index, and the short hand notation for gradients 
Bj = Vj-B). Here ai are the viscosity (Leslie) coefficients, A represents the symmetric part 
of the fluid velocity gradients [Aij = \{vij + Vj t i)] and the vector Ni = hi + ~ [n x curlv]^ 
is the change of the director with respect to the background fluid. Finally, there is a static 
(elastic) contribution due to the curvature of the director field 

o% = -Kn kyj n k)i , (2) 

given in the one constant approximation (Frank elastic constants K\ = K% = K3 = K). 

The director field is determined by the balance between the static molecular field h° = 
KV 2 h and the viscous molecular field hl i = (a^ — as)Ni + (a^ — a^jrijAij . The total molecular 
field has to be parallel to the director but h! can be neglected in the low Ericksen number 
regime (Er = avR/K <C 1, where v is a characteristic velocity and R the radius of the 
sphere). This condition is met in a typical thermotropic liquid crystal with K ~ 10~ n N; a ~ 
(5 — 10) x 10 _2 Pa • s in the case of v R 10~ 8 m 2 s _1 which allows speeds of millimetre per 
second for small colloid particles (R ~ 10/im). 

Considering low Reynolds number flow and using the equation of continuity we end up 
with seven equations 

a ihj = (3) 
Vi,i = (4) 
Kriijj = Xrii (5) 

for seven unknown variables (three for the velocity field v, three for the director n and the 
Lagrange multiplier A constraining n 2 = 1, and one for the pressure). The equation ([|) for 
the director is decoupled from the velocity due to the low Ericksen number approach and 
can be solved separately for the static problem, which then leaves only the hydrodynamic 
part of Eqs.(§|). 
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Once the velocity field v(r) and the pressure p(r) are obtained, the convenient way to 
determine the drag force is by calculating the total dissipation in the system 

F-v oc = J (tr' : A + h' ■ N) dV (6) 

where is the constant velocity of the fluid at infinity. 

III. NUMERICAL METHOD 

We followed the example of Heuer, Kneppe and Schneider || and used the artificial 
compressibility method (see for example ||) to solve the equations of motion. The idea of 
this method is that the system starts with an arbitrary start-up velocity and pressure field 
and relaxes in an artificial time towards its equilibrium, which is the solution we are looking 
for (d t p = 0, d t Vi = 0). The equations to solve are: 

v {ji = -c 2 d t p (7) 

where c is an arbitrary damping parameter, which should be chosen as large as possible to 
speed up the calculation (however, if c is too large the numerical scheme becomes unstable). 

Due to the linearity of Eqs.(|7]) it is necessary to solve them only for the two particular 
cases where flow and director are parallel and perpendicular to each other [n(oo) || v(oo) and 
n(oo) _L v(oo)]. The advantage of these solutions is the simple geometry. For an arbitrary 
angle between velocity and director they are just added together, i.e. the friction drag can 
be calculated by the resistance tensor M^- = M±5ij + (M\\ — M±)niTij, which determines the 
response of the drag force on the sphere to the flow around it: 

Fi = Mi^h^j = M ±Vooi + (My - M ± )( Voo ■ n)m. 

For an isotropic liquid the tensor is simply My = M8ij where the constant M is given 
by the Stokes friction M = —QtiRt]. The ratio Mj_/M\\ is a measure for the anisotropy 
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in the system since this gives the lift effect in the drag force. In the first case, assuming 
that the flow is along the z-axis and the director is parallel to it, the system is symmetric 
with respect to azimuthal rotations around the z-axis. When the velocity components are 
transformed to cylindrical components (v x ,v y ,v z — > v p ^v^v z ) the azimuthal velocity is 
zero everywhere in the system and the problem becomes two-dimensional. Furthermore it is 
favourable to use a spherical coordinate system with an inverse radius (£ = - = —, — = ; 



[v(oo), p(oo)] are included in the grid used for the calculations and the mesh size of the 
grid is smaller near the surface of the sphere, where most of the changes happen, and large 
far from the particle, where the values stay almost constant. It is sufficient to pursue the 
calculations in one quadrant only since the other three are given by symmetry. It is also 
evident that the radial velocity must be zero at both boundaries. The values for v z at 
9 — | can be computed as the inner grid points whereas the values at 9 = request special 
treatment since they contain the term cosec 9 and are therefore of the form "0/0". Since it 
was not possible to obtain them by an interpolation, we simplified the equation by taking 
the limit for 9 = [v z (£,Q) = lim^o Vz{£, 0)] analytically (application of L'Hopital's rule). 

In the second case, the director perpendicular to the velocity, there is no rotational 
symmetry and the calculations have to be done on a three dimensional grid. It is 
again favourable to use spherical coordinates with an inverse radius (see fig.|2|) for the 
reasons explained above, but this time the velocity components are kept Cartesian [i.e. 
v x (£, 0, Q),v y (£, 4>, 9),v z (£, 0, 9)]. Due to the symmetry it is sufficient to solve the equations 
in one octant. The conditions on the boundaries of this octant and the needed values are 
calculated as is shown in table 1. 

The constant number of grid points in the plane of the azimuthal angle direction 
(independent of 9) leads to a decrease in the mesh-size in real space while approaching the 
pole and finally yields a non- uniqueness for the pole (9 = 0) itself. A constant distance in 
real space would request fewer points (factor ~ 0.7) but it involves more calculations since 



9 = arctan 



z 
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the derivatives become more difficult. Therefore we chose the grid shown in fig.|2]. The values 
at the z-axis, which are non-unique, were then calculated for each <fi and set to their average 
over <p {v z {£, 0, 0) =< v„(g, (j), 0) >< t> for every f ). 



IV. THE DIRECTOR FIELD 

As described above, the director field can be taken as fixed during the calculations in 
the low Ericksen number regime. In order to study the influence of simplifying assumptions 
concerning the form of the field we performed the calculations with different director fields 
n(r). In the one constant approximation the director field is described by the minimum of 
the Frank free energy Fd 

F d = J (V • n) 2 + (Vn)W . (8) 

If we take into account that n is a unit vector and set n(oo) parallel to the z-axis, we can 
write the director components as 

n x = sin (3 sin 7 (9) 
n y = sin/3 cos 7 (10) 
n z = cos (3 (11) 

where (3 and 7 are angles dependent on the spatial coordinates. The director field in our 
problem is rotational symmetric with respect to the z-axis. We can therefore set 7 = 
arctanjz/x. Inserting this in Eq.(§) and minimising it, we are left with one equation for the 
polar angle /3(r): 

V 2 /J-^L=° (12) 
2r z sin V 

where (3 is a function of the radius r and the azimuthal angle 9. There are several possibilities 
to proceed with finding the static director field (see fig. 3): 
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If we neglect the boundary conditions on the surface (anchoring energy is zero) we get 
the (trivial) solution /3(r,6) = 0, i.e. the director is uniform in space, parallel to the 
to the z-axis. This was the approach chosen by the authors H in their analysis. 

Provided the anchoring on the surface is weak and therefore the angle (3 [the deviation 
from n(oo)] remains small, Eq. (|r2|) can be linearised and yields 



V 2 /?-^- = 0. (13) 
r z sin 6 



This equation, with the corresponding boundary conditions and the symmetry of the 
problem, can be easily solved and gives (3 = (j^rj R 3 sin2#/r 3 , where W is the an- 
choring energy and R the radius of the particle 0. 

If we assume strong anchoring on the surface of the sphere (WR/K ^> 1) the director 
field is forced to have a disclination loop (radius a) around the equator of the field. 
The direction of the molecules in the plane of the ring must be radial between ring 
and disclination, and parallel to the z-axis in this plane outside the disclination. Fur- 
thermore the perturbation of the field must decay as 1/r 3 far from the sphere. The 
simplest function which shows this behaviour is 

(3 = arctan g. (14) 

2 cos 20 + (*) 

There is an extensive discussion of the features and details of the director field in 
the strong anchoring regime |7[]. The conclusion reached [[!]] is that Eq.(|14D provides 
a very good approximation, well describing the far-field behaviour, the disclination 
ring vicinity, and even the weak anchoring case when the ring radius a is taken a — > 
WR 4 /AK. 



There is no analytical solution for the whole problem (satisfying Eq. (|T2|) and the 
boundary conditions). We therefore solved the equilibrium equation ( |i2| ) numerically 
with a method similar to the artificial compressibility method mentioned before. In 
this way we obtain the exact director field n(r) on every point of our grid. 
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V. RESULTS 



First we examine the influence of the different director fields [uniform, trial function f3 
and the exact numerical n(r)] on the drag force acting on the sphere, using the particular 
set of viscous coefficients of MBBA [ffl. As expected, the uniform field shows a much lower 
drag force for both principal configurations of n(oo) and v(oo) than the trial field and the 
exact field, see fig. 4 (for the velocity parallel to the director it is even smaller than for the 
isotropic drag force). On the other hand the anisotropy of the drag forces (the ratio of 
the two forces F±/F\\) becomes smaller for the realistic non- uniform field. This is due the 
following: the flow velocity around the sphere is the highest in the region of the equator 
plane perpendicular to the line of motion. If the liquid crystal is oriented along the same 
z-axis we get a very high gradient of the director in exactly the same region due to the 
effect of the disclination ring. In the other principal configuration, where the director n(oo) 
is perpendicular to the z-axis and, therefore, perpendicular to the velocity, the disclination 
with its high gradients is the same, but this time the loop is around a longitude of the 
sphere. It still increases the drag but in a much smaller region since the flow velocity at the 
stagnant poles is almost zero already. The anisotropy in the drag force for the three director 
fields yield: 

F j_ F* _i_ F* j_ 

-^H uniform = 1-69 "^H trial/3 = 1-50 —(exact = 1-50 (15) 

The results obtained with the trial field are surprisingly close to those of the exact 
director field. They obviously reflect the important features of the field which are mainly 
the disclination loop and the 1/r 3 decay of the deviation in the angle far from the particle, 
whereas the particular details near the particle and the disclination seem to be of minor 
importance. Therefore, the drag force is determined by the long-range effects. The difference 
between the drag force of the trial functions and the drag force obtained from the exact 
numerical solution is less than 1% and the difference in the force ratios is smaller than the 
accuracy of the calculations. In most practical cases it should be sufficient accuracy to 
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use these trial n(r) fields instead of a numerical solution of the governing equations. The 
uniform field, on the other hand, is not a very useful assumption since its resulting drag force 
differs from that for the exact n(r) distribution by up to 20%. The significantly higher ratio 
is a particular problem since it shows that the effect of boundary conditions on the surface 
cannot be simply modelled by a larger effective hydrodynamic radius. The dependence on 
the mutual orientation of rzoo and Voo is too strong. 

These results can be compared with experimental figures. The diffusion of particles 
in a nematic liquid crystals |§ is described by the diffusion tensor D which is indirectly 
proportional to the resistance tensor D = /cbT(M) -1 , where ks is the Boltzmann factor 
and T the temperature. Consequently, the tensor is of the same form as the mobility tensor 
Dij = Dj_5ij + (.Dm — D±)ni7ij with Dj_ = ksT/Mj_ and Dii = ksT/M\\. Recent experiments 
H showed for the the self diffusion constants of MBBA a ratio of D\\/D± ~ 1.5. Easier to 
determine experimentally is the anisotropy in the electric conductivity [|12j /x of a sample, 
which is related to the diffusion by 

D = — -/x 16 

ne z 

for n charge carriers of charge e per cm 3 . The conductance anisotropy was often measured 
H] and for MBBA it is usually equal to ~ 1-5. Both experiments are in excellent 

agreement with our result for Mj_/Mh = 1.50. 

The dependence of the drag force on the temperature is also of great interest in many ex- 
periments. The viscous coefficients a±, a 2 , a 3 , a 5 and a 6 , in the first approximation, depend 
linearly on the order parameter S in the region close to the nematic^isotropic transition 
temperature T ni . The order parameter S itself can be approximated by Haller's equation 




S=(l-^-Y (17) 
where 7 is determined experimentally for MBBA [O] to be 7 = 0.188. The viscous coeffi- 



cients scale, therefore: 
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cui — > «i * 5; «2 — ► <^2 * S*; 0:3 — ► 0:3 * S 1 ; 0:4 — > 04; 0:5 — * 0:5 * 5; ai6 — > ai6 * 5*. 

(18) 

The drag force shows for the perpendicular case more or less the same behaviour for 
all director fields: after a jump at the transition point, it increases while lowering the 
temperature, in the beginning rapidly, then slower and slower. On the other hand there is 
a qualitative change for the parallel drag force: while it jumps to a lower value and then 
decreases further for the uniform field, it shows a small change to a higher value at the 
transition point at which it stays almost constant independent of the temperature, both for 
the trial field of (3 and the exact field (see fig. 4). 

The boundary conditions are not absolutely rigid in many experiments due to a finite 



anchoring energy, in which case the approximative director field (|i~3D can be used. A typical 
colloidal particle of radius R = 10~ 5 m in a liquid crystal with an elastic constant of K ~ 
10~ n N has an anchoring energy of W ~ 10~ 5 . . . 1CT 7 J. That corresponds to a relevant 
dimensionless factor of WR/K = 0.1 ... 10. Our calculations showed a slow linear increase 
of the drag drag forces F± and Fu in the range of weak anchoring, WR/K = ... 4 due to 
the deviation of the director field from the uniform state. This results in gradients in the 
field which increase the dissipation and, therefore, the resistance of the particle to the flow. 
The effect on the force in the parallel case Fu is obviously stronger than in the perpendicular 
case F± which is reflected in a linear decrease of the anisotropy ratio F±/F\\ while enlarging 
the anchoring energy W. 

The authors of have also examined the case of charged particles, when the radial 
electric field near the surface forces the disclination loop to be pushed further away from 
the particle. An approximate expression for the loop radius a is then given by 

2 _ e a ( l 2 /in\ 

~ 32e 2 (5K + K 13 -2K 24 ) 1 ' 

where e a is the dielectric constant, q the charge of the sphere, and K13 and K24, elastic 
constants (assuming no immediate screening). One expects an increase in the drag forces 
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and a decrease in their ratio since the hydrodynamical effective cross section of large director 
gradients is far more increased for the parallel case than for the perpendicular one. For 
instance, taking the loop radius a ~ 2R the results for the MBBA-set of Leslie coefficients 
are: 



F ± = 1.75F U 



1.23 



F± 
Fa 



1.43 



(20) 



where the drag forces are given in units of MBBA in the isotropic phase F iso = —6nRr]v 
using rj = 0.5^4 as viscosity coefficient. As mentioned above, the anisotropy of the drag 
forces decreases with increasing the strength of the boundary conditions from F±/ Fh—1.69 
for the uniform field (anchoring energy W — 0), followed by a slow linear increase for weak 
anchoring (WR/K -C 1), and F±/F\\=1.5 for rigid anchoring to F±/F\\=1A3 for the case of 
the charged particle, which can be considered as "over-strong" anchoring. 

The molecular characteristics of the liquid crystal are inherent in the viscous coefficients. 
These coefficients depend, among other things, on the shape of the molecules which form 
the liquid crystal. This influence can be modelled by an affine transformation model [13" 
giving the viscous coefficients depending on the molecular aspect ratio £\\/£±: 



1 £\\ £ ± " 

ai = -2 ao {£;-j;. 



a 2 = ( 1 - 



(21) 



"3 = 



-± 



Q!4 = «4 



a 5 = —a 2 



«6 = «3 



The authors of |L3| determined the aspect ratio of MBBA to be £\\/£± = 5/2. We used this 
value to obtain the constant cto by comparing the largest coefficient a 2 with the experimental 
value and the isotropic coefficient was taken from MBBA directly. We calculated the cases 
of two particular configurations: a rod-like molecule with ratio £\\/£± = 7/2 and a disc-like 
system £\\/£± = 3/5. The results of these calculations are: 



F± 

1.78 F iso "77" 



rod-like molecule: Fx = 2.33F iso i 7 
disc-like molecule: F ± = 0.94 F iso F\\ = 1.52 F iso 
Note the inverted ratio of the drag forces for the disc-like molecules. 



F± 
Fi 



1.31 
0.62 



(22) 
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VI. CONCLUSION 



Considering the low Ericksen number regime (Er = avR/K <C 1) the director field can 
be taken as independent of the flow in the first approximation. We, therefore, took several 
static director fields (approximations for the field in the limit of strong and weak anchoring 
and the solution of the governing equations). The equations of motion were then solved 
numerically for the different fields n(r) using the viscous coefficients of MBBA. Due to the 
linearity of the Eqs.(|7]) it is only necessary to solve two limiting case, for director and velocity 
parallel and perpendicular to each other at infinity [n(oo) || v(oo) and fi(oo) _L v(oo)]. This 
yields the drag forces F\\ and F± which can be combined for the general case by the mobility 
tensor. 

The comparison of the drag forces for the different director fields showed that the discli- 
nation loop around the sphere, which is topologically necessary for a large anchoring energy 
of the molecules on the surface, does not only increase the forces itself but also decreases 
their ratio to F±/Fu = 1.50 compared to a uniform case (F±/Fu = 1.69). Trial director 
fields, constructed from the basic features of the director field (disclination ring and 1/r 3 
decay of the far field), showed to be a very good approximation. The difference between the 
values of F\\ and F± compared to the ones obtained for the exact field is less than 1%. 

The temperature dependence of the drag force showed an increase in the force F± for 
decreasing temperature and an almost constant value for the parallel force (F\\). An approx- 
imation for the director field for weak anchoring energy shows a linear decrease for lowering 
the anchoring energy in both particular drag forces as well as in their their ratio. 

The value for MBBA (F±/F» = 1.50), using the exact solution of the director field, is in 
good agreement with experimental results, measured by the static conductivity and the self 
diffusion of MBBA (both ~ 1.5). 

The disclination loop can be pushed away from the sphere, for instance, in the case of 
non-screened charges on the particle. This increases the particular forces compared with the 
uncharged case, where the loop is close to the surface, and yields for a loop radius of twice 
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the particle radius an even lower ratio of F±/F\\ = 1.43. 

The viscous coefficients of other materials can be approximated by an affine transforma- 
tion model, which uses the aspect ratio (£\\/£±) of the molecules as parameter. For a rod-like 
molecule (£\\/£±=7/2) we obtained an anisotropy in the drag force of F±/F\\ = 1.31 and for 
a disc-like molecule (i\\/£±=3/5) the ratio obtained was F±/F» = 0.62. The ratio smaller 
than one indicates that the lift force turns the particle away from the director whereas a 
ratio larger than one forces the particle in the direction of the director. 

It is especially interesting to examine the lift component of the drag force, i.e. the non- 
dissipative force acting perpendicular to the line of particle motion. It resembles magnetic 
forces and leads to physical phenomena, similar to the Hall effect. In a long cell, the ratio 
between the cross voltage U* and the applied voltage U is determined by the anisotropy 



ratio V = fi\\/fi± |T2 



— = b Sin29 (23) 

U a(V + l)/(V- 1) -cos2# 1 ' 

where a is the width of the sample (in direction of the applied voltage) and b the thickness 
of the sample (in direction of the cross voltage). The conductivity is determined by the 
movement of the charge carriers and, therefore, inversely proportional to the resistance 
which yields for the anisotropy ratios V = M±/M\i = F±/Fu. 

Since it is difficult to produce samples of liquid crystals without disclinations, which are 
sufficiently large to perform measurements of moving particles, the effect of the lift force (the 
component of the drag force which is perpendicular to the driving force) could be observed 
in a long cylinder. If the boundaries force the liquid crystal to be perpendicular to the 
walls of the cylinder, the director field will "escape to the third dimension" , i.e. it will turn 
round to be parallel to the long axis of the cylinder while approaching its centre since this is 
energetically much more favourable than a disclination line. If the distances with the same, 
well defined curvature along the z-axis are large enough, little spheres, which are dropped 
in the sample, should show a certain measurable displacement during their way down in a 
gravitational field. 
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Further possibilities are the usage of electric and magnetic fields. Moving particles can 
be guided by changing the director orientation in the sample to direct them to a certain 
destination in the sample. This enables the guiding of uncharged and unpolarizable particles 
with electric or magnetic fields 
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Fig.l: If a sphere falls down in a gravitational field and the director field is not parallel 
to the force, their is besides the friction F p anti-parallel to the gravitational force also a 
component Ft, perpendicular to it, the so-called lift force. Note also the disclination loop 
around the particle due to the boundary conditions, indicated by the dots. 

Fig. 2: One layer of the matrix H = \ =const.) and its transformation to real space. 
The lines of constant angle transform to longitudes and constant 9 to latitudes. Note 
the decreasing distance between two longitudes while approaching the pole. This yields a 
non-unique point at 9 — 0. 

Fig. 3: The director fields for a) no anchoring at the surface of the sphere b) weak 
anchoring c) arctan field d) numerical solution fulfilling the boundary conditions. The black 
dots show the position of the disclination loop, the grey dots show where they would be 
(having this director field without the sphere). 

Fig. 4: The drag force on the sphere for three different director fields (uniform, trial and 
real field) depending on the effective temperature T/T ni . The upper lines are for the case 
of director and velocity parallel , the lower ones for director and velocity perpendicular to 
each other. 

Table 1: The treatment of the boundaries around the first octant. The notation be refers 
to the fact that this value is fixed by the boundary conditions of the problem, symm refers to 
given by symmetry, the values given by cal can be calculated with the untreated equations 
of motion, se/inter the value was determined by an analytically simplified equation for the 
uniform field and by an interpolation for the non-uniform field, finally extr means the value 
was calculated by an extrapolation. 



23 



